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ABSTRACT 

We present the discovery of a large-scale structure of emission-line galaxies at 
redshift z = 4.86 behind a massive cluster of galaxies, A1689. Previous spectroscopic 
observations of a galaxy, A1689-7.1 at z = 4.87, near this structure, revealed a possible 
overdense region of inter-galactic medium (IGM) around the galaxy, which extends at 
least ~ 80 comoving Mpc along the line of sight. In order to investigate whether this 
z b IGM overdense region contains a galaxy overdensity, we undertook narrow- 
and broad-band imaging observations around A1689-7.1 with Subaru/Suprime-Cam. 
We detected 51 candidate Lya emitters (LAEs) at redshift z — 4.86 ± 0.03 in the 
32 X 24 arcmin^ field of view. After correction for lensing by the foreground cluster, 
we found a large-scale 20 x 60 comoving Mpc) overdense region of galaxies around 
A1689-7.1 in the source plane at z = 4.86. The densest peak in this region has an 
overdensity of 5 4, suggesting that this structure is probably a good candidate for 
a protocluster which may evolve into a massive cluster of galaxies in the present-day 
Universe. A1689-7.1 is located at the edge of this region, where the local galaxy density 
is ~ 1.6 times the mean density and is close to the density contrast in the IGM along 
the line of sight to A1689-7.1 estimated from the optical depth. The overdensities of 
galaxies we have found may suggest that at least some parts of the IGM overdense 
region have already started to form galaxies and moreover they relate to the formation 
of a protocluster. Although we lack information on the three dimensional distributions 
of both IGM and galaxy overdense regions, the similarity of the scales of both regions 
may suggest that the two are parts of a single large-scale structure, which would be 
an large edge-on sheet along the line of sight with a size of ^ 20 x 60 x 80 comoving 
Mpc. 

Key words: cosmology: observations ~ galaxies: formation ~ galaxies: evolution - 
galaxies: individual: A1689-7.1. 



1 INTRODUCTION 

Galaxies are believed to interact frequently w ith the sur- 
roun d ing intergalactic medium (IGM) ( e.g., iKeres et al.| 
l2005l : iKobavashi. Springel. fc White! l2007l ). There are sev- 

* Based on data collected at Subaru Telescope and in part ob- 
tained from data archive at Astronomical Data Analysis Center, 
which are operated by the National Astronomical Observatory of 
Japan. 
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eral observational studies of these interactions using corre- 
lation between_gal^jdes and IGM a b sorption- lines a t z < 3 
(e.g., lAdelberger et al.l 12001 l2005l : ICrighton et al] |2009| ). 
While the studies of IGM absorption-lines usually need a 
bright backgrou nd source, such as a QSO o r a gamma- 
ray burst (e.g., lRauchlll998l : IChen et alj|2005l ). iFrve et all 
(2008) claimed the existence of an overdense region of the 
IGM at 2; ~ 5 using a galaxy, A1689-7.1 a.t z — 4.866, which 
is lensed by the foreground massive galaxy cluster A1689 at 
z = 0.187. 
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Lya Redshift 
4.70 4.75 4.80 4.85 4.90 4.95 5.00 



Table 1. Summary of Observations 
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Figure 1. Transmission curves of NB704 and NB711 filters. The 
horizontal bar and vertical thick bar indicate the Lya wavelength 
of the HI overdense region and the redshift of the l ensed galaxy, 
A1689-7.1, 2 = 4.866, respectively jFrve et al.ll2008l'l . 



In the spectrum of A1689-7.1, the HI flux transmission 
at 6920-7100 A (Lya HI absorption- lines at z = 4.69-4.84) 
is only ~ 30% of the average value at the same redshift de- 
rived from z > 5 QSO spectra, suggesting that this over- 
density of inter-galactic H I gas clouds around A1689-7.1 
extends at least ~80 comoving Mpc along the line of sight. 
Although individual H I absorption-lines could not be re- 
solved at the resolution of their spectrum, they claimed that 
this structure is likely to be made by overlapping of many 
intergalactic H I gas clouds rather than a few overlapping of 
very high H I column density clouds, because of absence of 
metal absorption-lines from these clouds on the continuum 
long-ward of the Lya emission. 

Such a large-scale IGM overdense region would provide 
us with an unique opportunity to study the formation of 
galaxies within baryonic large-scale structure at z ~ 5. Al- 
though there have been several high-z emission-line galaxy 
searches targeting QSO absorption- line systems, the targets 
for most of these searches were absorption- line systems by 
very high H I column density clouds, such as damped Lya 
absorbers rather than large-scal e overdense regions of low H 
I col u mn density c l ouds (e.g.. iFvnbo. MoUer. fc Thomsen 
I2OOII: iFrancis et al.l liooH : iPalunas et all |2004| ). iFrve etal 



1 20081 ) claimed that this IGM structure may be an exam- 
ple of the earliest phase of baryonic large-scale structure 
that subsequently evolve into the filamentary cosmic web 
of galaxies (including clusters of galaxies) observed in the 
present-day Universe. How can we test the possibility? Does 
the IGM overdense region already contain galaxies and even 
galaxy overdensities? If the overdensity of galaxies exists, 
what is the size, density, and morphology of the structure? 
In order to answer these questions, we undertook a wide- 
field, deep Lya imaging survey for emission-line galaxies at 
z ~ 4.8 in the A1689 field. 

In this letter, we use AB magnitudes and adopt cos- 
mological parameters, JIm ~ 0.3, Qa ~ 0.7 and Ho — 70 
km s~^ Mpc~^. In this cosmology, the Universe at z = 4.86 
is 1.2 Gyr old and 1.0 arcmi n corresponds to a comoving 
length of 2.3 Mpc at z = 4.86 l|Wrightll200^ ). 
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"The central wavelength and FWHM of the filters. 
'The 5cr limiting magnitudes within 1.6 arcsec diameter aper- 
tures. 

'^Significant ghosts appeared in NB70A images. 



2 OBSERVATIONS AND DATA REDUCTION 

We obtained narrow- and broad-band imaging centred at 
{a,S) = 13:11:29.5, -01:20:17 (J2000.0) on 2009 M ay 23 - 
24 (UT) using Suprime-Cam (Mi vazaki et al.ll2002l ) on the 
8.2-m Subaru Telescope (lyc et al.i 20041 ) as an open use nor- 
mal programme (S09A-104). We also used archival V'-band 
data of the same field taken on 2001 Apr 24 (UT). De- 
tails of the observations are listed in Table 1. Our original 
plan was to use tw o narrow-band filters, N B70A and NB711 
(|Ouchi et al.ll2003l : [Shimasaku et al. 2oM [200i . see Fig. 1). 
We prioritized NB704 imaging because NB7QA covers the 
redshift range for Lya at z ~ 4.76 — 4.84 (46 comoving Mpc), 
which perfectly matches the redshift range of the IGM over- 
dense region at z ~ 4.69 — 4.84. Unfortunately, the filter 
holder of NB704: had problems and significant ghosts ap- 
peared in our scientific images We immediately changed 
filter from iVB704 to iVB711. Although the redshift cov- 
erage of iVB711 (z ~ 4.83 - 4.89, 33 comoving Mpc) has 
a smaller overlap with the redshift range of the IGM over- 
dense region than NB7Q4:, NB711 imaging is still useful to 
search for emission-line galaxies around A1689-7.L 

The raw d ata were reduced with SDFRED (|Yagi et al.l 
l2002l : lOuchi et a l. 2004) and iraf. We flat fielded using the 
median sky image after masking objects. We did background 
sky subtraction adopting the mesh size parameter of 64 pix- 
els (13 arcsec) before combining the images. We calibrated 
the astrometry of the i mages u sing the 2MASS All-Sky Cat- 
alog of Point Sources (|Gutri et al. 2003,). Photometric cali- 
bration was obtained from the photometric and spectro pho- 
tome t ric standard st a rs, SA104, SAllO-2 32, and Hz44 (|Okel 
I1990I : iLandoltl 1 19921 : ISmith et all I2OO2I ). The magnitudes 
were corrected for Galactic extinction of E{B — V) — 0.03 
mag ((Schlegel. Finkbeiner. fc Davislll998l ). The variation of 
the extinction in this field is small (peak to peak, ±0.006 
mag) and thus it does not affect our results. 

The combined images were aligned and smoothed with 
Gaussian kernels to match their seeing to a FWHM of 0.8 
arcsec. We made a Ri image [Ri — {R + i')/2] for the contin- 
uum with an effective central wavelength similar to NB711. 
The total size of the field analyzed here is 32.1 x 23.8 arcmin^ 
after removal of low S/N regions near the edges of the im- 



^ The filter holder problems were fixed by the staff of the Subaru 
Telescope immediately after the run. 



Galaxies and HI gas at z = 5 3 



4 



1^ 
m 




22 23 24 25 26 

NB71 1 (ABmag) 



Figure 2. Colour-magnitude plot of NB711 vs Ri - NB711 
for NB711 detected sample. The solid line indicates colour and 
magnitude limits {Ri - NB711 = 0.8 & iVB711 = 25.7). The 
red circles indicate 51 candidate LAEs. The green square indi- 
cates A1689-7.1. The blue curve indicates the 3a uncertainty of 
Ri — NB711 colour based on photometric errors of both Ri and 
A'^B711 for objects with constant spectra. The emission-line 
objects that are not circled are considered to be low-z interlop- 
ers from their broad- band colours (see text). All magnitudes and 
colours are measured with 1.6 arcscc diameter apertures. 

ages. We also masked out the halos of the bright stars. The 
resultant total effective area is 723 arcmin^. 

Object detection and ph otometry were peri:ormed using 
SExTRACTOR version 2.5.0 iBertin fc Arnoutd Il996l ). The 
object detections were made on the NB711 image using a 
Gaussian detection kernel with FWHM of 0.8 arcsec. We 
detected 173,864 objects that had 5 connected pixels above 
1.5a of the sky background. The magnitudes and colours are 
measured in 1.6 arcsec diameter apertures. The magnitude 
limits of the observations are reported in Table 1. 



3 RESULTS 

Fig. 2 shows a colour-magnitude plot for the A'^_B711 de- 
tected sample. We selected 51 candidates as Lya emitters 
(LAEs) down to iVB711 = 25.7 (5ct limit) with the follow- 
ing criteria; 

(1) Ri - iVB711 > 0.8 (EW^bs > 85 A), 

(2) R-i' > 0.5, 

(3) i' -NB7n > 0, 

(4) V - i' > 1.7 or V > 27.3 (2a) , 

(5) E > 3, 

where the E is the ratio between the A'^_B711 excess and the 
uncertainty of Ri — NB711 colour based on photometric er- 
rors of both Ri and A'^_B711 for objects with constant f,j 
sp ectra. The cr i teria, (1), (2), and (3) are the sam e as used 
bv lOuchi etal] (|2003l ) and lShimasaku et all l|2003D to select 
2 = 4.86 LAEs with the same filter sets. The criterion (4) is 
apphed to exclude the contamination from low-z emitters, 
such as [OII]A3727 emitters &t z = 0.91, [OIII]A5007 emit- 
ters at z ^ 0.42, and Ha emitters at z = 0.09. The V — i' 
colour limit in (4) wa s determined us ing the expected colour 
of z = 4.9 galaxies bv lSongailal l|2004 ). The criterion (5) was 
applied to exclude the contamination from non-emitters due 
to photometric errors. The properties of the 51 LAEs are 



listed in Table Al and displayed in Fig. 3. We note that 
A1689-7.1 was not selected as a LAE, because it failed to 
satisfy the criteria, (1), (3), and (5). 

We corrected the measured magnitudes of LAEs to 
intrinsic (un lensed) magnitudes us ing the mass model of 
A1689 from |l imousm et all (|2007h which combines both 
stron g- and weak-lensing constraints (see also I Richard et al] 
I2OO9'). While the ma gnification facto r of A1689-7.1 equates 
to 2.5 magnitudes (jFrve et al.l 120081 ). for our LAE sample, 
individual magnification factors range from to 0.9 mag- 
nitudes. We found that the intrinsic A'^i3711 magnitudes 
{NB711int) of six LAEs are fainter than the iV_B711 mag- 
nitude limit (A'^i3711 = 25.7). Lensing also reduces our sur- 
veyed area by the magnification factor, and we estimated a 
surface of 665 arcmin'^ in the source plane, 8% smaller than 
the observed field of view. We also estimated from the mass 
model the angular displacement of each candidate back to 
the source plane at 2 = 4.86, typically 10 to 60 arcsec to- 
wards the centre of the cluster (the displacements from the 
measured positions, see Fig. 4) and perform all further cal- 
culations in the source plane. 

We tested whether the entire field of view of the A1689 
image (in the source plane) exhibits an overdensity by com- 
paring to the A^_B711 surveys in Subaru Deep Field (SDF) 
llOuchi et al ] |2003l : IShimasaku et al ]|2003l ). In the SDF, they 
detected 43 LAEs down to NB711 = 25.5 in 25x45 arcmin^. 
The number density of LAEs in the SDF is 0.04 ± 0.01 
arcmin"'^. The uncertainty is a la Poisson error. Applying 
the same (intrinsic) magnitude limit to the A1689 field data, 
we detected 33 LAEs in 665 arcmin^. The number density 
of LAEs in the A1689 field is 0.05 ± 0.01 arcmin"^ Thus 
the number density of LAEs in the A1689 field is similar to 
that in the SDF at least for a sub-sample of slightly bright 
LAEs (L(Lya)> 1.7 x 10**^ erg s"^). Thus the A1689 field 
does not show evidence for overdensity in the entire field of 
view if we use the SDF as a blank field. 

Fig. 4 shows the spatial distribution of the 51 candidate 
LAEs in the A1689 field and the smoothed density map in 
the source plane. To avoid effects from the spatial distribu- 
tion of the lensing magnification, we made the density map 
using intrinsic magnitude limited sample of 45 LAEs with 
NB711i„t < 25.7 rather than the full 51 LAE sample. We 
used a top-hat smoothing kernel of 8 comoving Mpc radius , 
which is the same as that used in IShimasaku et al.l (|2003h . 
We defined high-density regions as the regions with the over- 
density 5=(n — n)/n>0. There is a large-scale overdense 
region with an extent of ~ 8 x 25 arcmin^ (~ 20 x 60 comov- 
ing Mpc) surrounding A1689-7.1. While the densest peak in 
this overdense region has an overdensity of 5 = 3.9, with a 
significance level of 4.3 a, A1689-7.1 is located at the edge 
of this region, where the overdensity is 5 = 0.6. 

We tested the reliability of the structure in the fol- 
lowing way: (a) We slightly changed the criteria for LAEs 
(i.e. ±0.1 mag), but the properties of the structure do not 
change significantly. We conclude that the structure is sta- 
ble against small changes in our sample definitions, (b) We 
divided the field of view into eastern and western parts, and 
selected low-z emitters from NB711 detected sources down 
to A''i3711 — 25.7, satisfying criteria (1), (3), and (5), but 
not satisfying (2) or (4) . The numbers of low-2 emitters in 
the eastern and western parts are 42 ± 6 and 55 ± 7, re- 
spectively. The uncertainties are la Poisson errors. These 
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Figure 3. True colour images {R for blue, Af_B711 for green, i' for red) of LAEs. The size of eacli imago is 8 X 8 arcsec'^ (51 X 51 kpc^ 
at 2 = 4.86 in absence of lensing). The LAEs are ranked in terms of intrinsic (unlensed) NB711 magnitudes. 




30.0 13:12:00.0 30.0 11:00.0 10:30.0 

RA (J2000) 

Figure 4. The sky distribution and smoothed density map (in the source plane) of candidate LAEs at 2 ~ 4.86 in the A1689 field. The 
gray scale image is the NB711 image. The thick horizontal bar indicates the angular scale of 20 comoving Mpc (3.4 Mpc in physical 
scale) at z = 4.86. The red and blue circles indicate candidate LAEs with intrinsic (unlensed) NB711 magnitudes, N BTllint ^ 25.7 
and N B7lli„t > 25.7, respectively. The sizes of the circles are proportional to measured emission-line fluxes. The black square indicates 
A1689-7.1. The lines from these points indicate the positions in the source plane at 2 = 4.86. The surface density map of LAEs (in the 
source plane) is generated with a top-hat smoothing kernel with a radius of 8 comoving Mpc. The black contours indicate deviations 
from the mean density of this field, 5 = {n — n)/n = 0, 1, 2, and 3, equivalent to significance levels of 0, 1.1, 2.2, and 3.3 a. The green 
contours indicate lensing magnifications of 0.01, 0.02, 0.05, 0.1, 0.2, 0.5, 1.0, 2.0, and 5.0 mag. 



numbers agree to within 30%, showing that there is no clear 
systematic eflFect to select NBUl excess objects, (c) We also 
selected red objects down to i' = 25.7, satisfying criteria (2) 
and (4). The numbers of red objects in the eastern and west- 
ern parts are 1403 ± 37 and 1097 ± 33, respectively. These 
numbers also agree to within 30%, showing that there is no 



clear systematic effect to select red objects with R — i' and 
V — i' colours. In contrast to low-z emitters and red objects, 
the numbers of LAEs in the eastern and western parts are 
8 ± 3 and 43 ± 7, respectivley. Thus systematic errors of 
colours in the field of view should not affect the large-scale 
structure of LAEs significantly. 
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4 DISCUSSION AND CONCLUSIONS 

What does the large-scale structure of galaxies (i.e., LAEs) 
we have found evolve into? This large-scale structure of 
galaxies is a good candidate for a protocluster. Numerical 
simulations of structure formation suggested that most (> 
80%) of structures with a size of ~ 20—40 comoving Mpc and 
with a peak overdensity of5j53atz~5 — 6 will evolve into 
rich clusters (M ^ lO^^/i~^M0) i n the present-day Un iverse 
|Suwa. Habe. fc Yoshika"TOll2006l : [o^erzier et al.ll2009l ). This 
structure has a size of ~ 20 x 60 comoving Mpc and the 
densest region has an overdensity of 5 ~ 4. Hence we con- 
clude that this structure is probably a good candidate for a 
protocluster. 

How does the galaxy overdensity relate to the IGM over- 
density in this structure? If we assume the optical depth in 
the IGM is proportional to the square of loc al gas densit y 
over the mean gas density [r oc (p/p)^] (e.g.. iRauchlllQQSh . 
then the gas density in the IGM overdense region would be 
estimated to be ~ 1.3 times th e mean gas densi ty at z 5 
by using the optical depth from lFrve et al.l (|2008f ). The local 
galaxy density at the position of A1689-7.1 is ~ 1.6 times 
the mean density, and hence it is close to that of the IGM. 
However, our current observations may have only probed 
the back-part of the IGM overdense region, and thus it is 
possible that the galaxy density will be higher in the fore- 
ground along the line of sight to A1689-7.1. Although there 
is a small overlap of the observed redshift ranges of the IGM 
and galaxy structures {Az ~ 0.01), the current results sug- 
gest that at least some parts of the IGM overdense region 
may have already started to form galaxies and moreover they 
are probably related to the formation of a protocluster. 

Is it possible that these galaxy and IGM overdense 
regions are parts of a single large-scale structure? Unfor- 
tunately, at present, we lack information on the three di- 
mensional distributions of both IGM and galaxy structures. 
Without spectroscopy, it is unclear whether the structure is 
concentrated in a particular redshift or smoothly extends to 
the whole redshift range {z ~ 4.83 — 4.89). Similarly, it is still 
unclear whether or not the IGM overdense region extends 
in a transverse direction or to higher redshift than that of 
A1689-7.1. Nevertheless, the similarity of the scales of both 
structures (a few times ten comoving Mpc) may suggest that 
the two are parts of a single structure. If the IGM and galaxy 
overdense regions are parts of a single structure, this would 
be an large edge-on sheet along the line of sight with a size 
of ~ 20 X 60 X 80 comoving Mpc. The future observations of 
LAEs at z ~ 4.79 with iVB704 will be a direct test of this 
possibility. 
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" The properties are based on photometry with 1.6 arcsec diameter apertures. Magnitudes below la are repla<;ed to la 

upper limits. 

^ The LAEs arc ranked in terms of intrinsic (unlensed) ArB711 magnitude, N Blllinf 
The measured positions of LAEs. 

Offsets (i5RA and 5Dec in arcsec) of LAEs on the source plane. The positive values show offsets to East and North 
directions. 



